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We investigate the distribution of occupation times for a particle undergoing a random walk 
among random energy traps and in the presence of a deterministic potential field U dct (x). When 
the distribution of energy traps is exponential with a width T g we find that the occupation time 
statistics behaves according to (i) the canonical Boltzmann theory when T > T g , (ii) while for 
T < T g they are distributed according to the Lamperti distribution with the asymmetry of the 
distribution determined by the Boltzmann factor exp(— U dct (x)/T g ) with T g and not T being the 
effective temperature. We explain how our results describe occupation times in other systems with 
quenched disorder, when the underlying partition function of the problem is a random variable 
distributed according to Levy statistics. 
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Consider a one dimensional Brownian motion, in a 
binding deterministic potential field U dct (x). Following 
the trajectory of an individual particle one may deter- 
mine the occupation time t , which is the time the 
particle spends in the domain x\ < x < x-i- In the 
limit of long measurement times t the occupation fraction 
p = t° cc /t is given by Boltzmann statistics, assuming 
that the dynamics is ergodic 
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where Z = exp[— C/ dct (a;) /T]dx is the normalizing 
partition function and T is the temperature. We see that 
the occupation fraction p does not depend on the details 
of the dynamics, for example on the diffusion constant 
or the initial conditions, which is of-course the strength 
and the generality of statistical mechanics. Majumdar 
and Comtet showed that when the potential energy is 
random the occupation fraction may exhibit large fluctu- 
ations from one sample of disorder to another, in particu- 
lar in the Sinai model was investigated. Generally, for 
disordered systems an important question is what is the 
distribution of the occupation fraction and its relation to 
the underlying disorder. 

Here we investigate occupation time statistics using 
the well known trap model 0, 0] ■ Briefly the model intro- 
duced in the seventies, describes a random walk among 
traps with a random depths of energy traps, with a den- 
sity of states which is exponentially distributed with a 
width T g (see details below). Such density of states leads 
to anomalous diffusion 0, [a] , and aging 0, HI, @, 01 when 
T < T g . The model, in fact a family of models, was used 
to describe dynamics of many systems: transport of elec- 
trons in amorphous materials [a, 0, [lfj| , single molecule 
pulling experiments [llj], rheology of soft matter [l^] e.g. 
emulsions, relaxation in glasses 0, 0, [l]| and green flu- 
orescent protein dynamics 14| to name a few. Recently 
the aging scenario of the trap model 0,03], was justified 
rigorously using the random energy model as the starting 



point [151 ]. 

Bouchuad Q introduced the concept of weak ergodic- 
ity breaking, in the context of the trap model, which 
implies the breakdown of Boltzmann's statistics when 
T < T g . The goal of this paper is to define and in- 
vestigate the rather strong deviations from Boltzmann's 
statistics in the quenched trap model. In particular we 
investigate the distribution of the occupation times for 
the quenched trap model quantifying the deviations from 
the standard canonical theory Eq. (JTJ) . At the end of the 
paper we will show how the main features of our theory 
can describe occupation times in several other models of 
quenched disorder, proving some generality of our results 
beyond the trap model. 

Quenched Trap Model We consider a particle under- 
going a one dimensional random walk on a quenched 
random energy landscape on a lattice. Lattice points 
are on x — 0, a, 2a, • • • , L where a is the lattice spacing. 
On each lattice point a random energy E x is assigned, 
which is minus the energy of the particle on site x, so 
E x > is the depth of a trap on site x. The traps ener- 
gies {E x } are independent identically distributed random 
variables, with a common probability density function 
(PDF) p(E) = (1/T g )exp(-E/Tg). Due to an interac- 
tion with a heat bath the particle may escape site x and 
jump to one of its nearest neighbors. The average time 
it takes the particle to escape from site x is given by Ar- 
rehhius law t x — exp(E x /T). Notice that small changes 
in E x leads to an exponential shift in t x . In particular it 
is easy to show that the PDF of the waiting times is 
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so when T <T g the average waiting time diverges. 

An additional bias is applied to the system. For ex- 
ample in transport processes this is a driving field, e.g. 
an external electric field (and see more details below). 
Let q x (1 — q x ) be the probability of jumping left (right) 
from site x respectively. The master equation for the 
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population on site x, P x is 
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For non biased random walks q x — 1/2 while for uni- 
formly biased random walks q x =/= 1/2 is a constant. The 
boundary conditions are reflecting q = and = 1 
though our main results are valid also for periodic bound- 
ary conditions. 

The local bias q x is controlled by a deterministic po- 
tential field U^ ot , which in some cases is controlled by 
the experimentalist. It is usually assumed that detailed 
balance conditions holds so that the dynamics of the pop- 
ulations reaches thermal equilibrium described by Boltz- 
mann's canonical ensemble. For the trap model this well 
known condition leads to 
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For example if a constant driving force field T acts on 
the system q x = 1/[1 + exp(^ r a/T)] [6|]. 

We consider a single realization of disorder in the ther- 
modynamic limit where the measurement time t — > oo 
before the system size is made large. One can show 
that the equilibrium of populations is described by Boltz- 
mann statistics, which is not surprising since we used the 
detailed balance condition. The total time the particle 
spends in the domain x% < x < L is the occupation time 
t° cc . This domain is called the observation domain. As- 
suming that the process is ergodic we have the occupation 
fraction for a single disordered system 
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is the partition function of the part of the system under 
observation and Z NO = £x=o° exp[-(L^ ot - E x )/T] is 
the partition function of the rest of the system. The oc- 
cupation fraction is a random variable which varies from 
one system to the other, the goal of this manuscript is to 
calculate its distribution. However first four comments 
are in place, (i) If we have only a single realization of 
disorder the occupation fraction is given by Boltzmann 
statistics, the question then is whether the occupation 
fraction a self averaging quantity. Namely is it repro- 
ducible in a second experiment when a different realiza- 
tion of disorder is investigated. More generally we have in 
mind the situation where one investigates many realiza- 
tions of disorder, for each the occupation fraction is ran- 
dom and hence one may construct its distribution, this 



case corresponds in principle to single molecule experi- 
ments where one may track independently a large number 
of individual molecules each one interacting with a unique 
random environment [l6| . (ii) We have assumed that for 
a single disordered system the dynamics is ergodic. This 
is so since we are considering the thermodynamic limit 
where the measurement time t — > oo first while main- 
taining a finite size of the system. For a finite system 
we always have a finite E max which is the maximum of 
the random energies {E x }, and hence there is always a 
long time t which is much larger than exp(_E max /T) after 
which the process is behaving according to the ergodic 
principle, (iii) The occupation time in Eq. ([5]) describes 
rather generally the occupation time of a particle in a 
random energy landscape and is not unique to the spe- 
cific dynamics of the quenched trap model. For example 
we could add random barriers to the dynamics of the 
model which would not alter its equilibrium, (iv) Pre- 
vious work [T3| considered the occupation times of the 
continuous time random walk (CTRW) model (annealed 
model), unlike the quenched trap model in the CTRW 
model ergodicity is broken and the system is not spa- 
tially disordered. 

From Eq. ([5]) we see that the distribution of the occu- 
pation fraction p is obtained in principle from the distri- 
butions of two independent random partition functions 
Z° and Z N °. Let G z o{z) and G z no(z) be the PDFs of 
Z° and Z NO respectively. Then the PDF of the occupa- 
tion fraction f(p) is found using Eq. ([5]) 
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We now consider the problem of finding G z o(z). 

If the deterministic part of the field U x ct is a constant, 
Eq. |6]) shows that Z° is a sum of independent identically 
distributed random variables, and then Gauss-Levy limit 
theorems apply. In contrast, when U x et is not a constant 
then we are dealing with the problem of summation of 
non-identically distributed random variables and hence 
in what follows we modify the familiar limit theorems for 
the case under investigation. 

Let n be the number of lattice points in the interval 
[a;i,i]. We consider the scaled random variable Z° = 
Z°/n 1 / Q with a = T/T g and T < T g . The Laplace 
z — > u transform of the PDF of Z° is found using Eq. 
© and p(E) 
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where ip(u) = exp(— ur)^(r)dr. We now consider 
the limit of large n. We use the small u expansion 
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where A — a\T(—a)\ and from Eqs. (|8I9[) we find 
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In the continuum limit of a — > n — > oo and L — x\ — 
an remaining finite we may replace the summation with 
integration and find the stretched exponential 
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where U dct (x) is the deterministic field in the continuum 
limit. The inverse Laplace transform of Eq. (fTTj) is the 
one sided Levy stable law. 

A similar calculation is made for Z N °. We invert the 
Laplace transform Eq. (fTTj) . switch back to the original 
variable Z° instead of the scaled one Z° , and find using 
Eq. 
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In Eq. (fT2"|) l a (z) is the one sided Levy stable PDF whose 
Laplace pair is l a (u) = exp(— u a ). Pb(T cj ) is Boltzmann's 
probability of finding the particle in the observation do- 
main calculated using the deterministic field with a tem- 
perature Tg 
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Solving the integral Eq. (fT2"|) we find the Lamperti [15 
PDF 
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Eq. p^l[T5l) are the main results of this manuscript, soon 
to be discussed in detail, which are valid in the glassy 
phase T <T g . 

For T > T g and in the same limit we have the usual 
canonical behavior 
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Eq. (I16| shows that when T > T g the disorder plays no 
role, indicating the reproducibility of Boltzmann's statis- 
tics Eq. |JJ, when the disorder is weak. 
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FIG. 1: The PDF of the occupation fraction for the de- 
terministic field U(x) — Tx. (a) For T = 3T 9 we find a 
delta function centered on the value given by Boltzmann's 
statistics. When T/T g — 0.7 [panel (b)] the non trivial dis- 
tribution of the occupation fraction has three peaks while for 
T/T g = 0.3 [panel (c)] the distribution is bi-modal. The + 
are simulations and the curve is the theoretical prediction Eq. 

l,T g = l,o = 10~ 5 and 



!fl5[) without fitting. We used T = 
the observation domain < x < 1. 



We now discuss the behaviors found in our main Eqs. 
(fT51 [T5)) . These equations give the distribution of the oc- 
cupation times, which is the generalization of the usual 
Boltzmann law Eq. (fTB"]) . The parameter TZ is called the 
asymmetry parameter, and if TZ = 1, f(p) is symmetric, 
for example when T/T g = 1/2 and TZ = 1 we get the 
arcsine PDF. The asymmetry parameter TZ is calculated 
by the usual type in integral over the Boltzmann factor, 
however now the temperature T g is the relevant tempera- 
ture not T [see Eq. (IT3l ] . Roughly speaking there are two 
sources of fluctuations: the disorder characterized by T g , 
and the temperature T. Hence when T <T g the relevant 
temperature is the "temperature of the disorder" that is 
T g . For example using Eq. (|13I15I16[) the average occu- 
pation fraction has the following surprising discontinuous 
behavior, 
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The average occupation fraction freezes in the colder 
glassy phase of T < T g in the sense that it does not de- 
pend on the temperature T, for any type of deterministic 
binding field. 

In Fig. [1] we demonstrate our results comparing our 
theory with numerical simulations on a lattice. We 
consider the situation where the deterministic field is 
U det (x) — Tx and < x, and the observation domain is 



< x < T g jT. In Fig. E^a) with T > T g we see that the 
distribution of occupation fraction is very narrow with 
p = P B (T = 3T g ) = 1 - e" 1 / 3 indicating that the dis- 
order is not important. In contrast when T < T g the 
behavior of the occupation fraction changes dramatically 
and p is non self averaging and random. 

Eq. {HjJ shows that when T/T g < 1 the PDF of occu- 
pation fraction is essentially composed of two delta func- 
tions centered on p — 1 and p = 0. Namely for some 
samples of disorder the particle is within the observa- 
tion zone during all the observation time t (p — 1) and 
in other samples the particle is never in the observation 
zone (p = 0). This behavior is easy to understand when 
T — > the minimum of the random potential energy is 
the most populated, and this minimum can be found ei- 
ther in the observation zone or out of it. As shown in Fig. 
[T](c), for small but finite T we have a non-trivial bi- modal 
U shape of the PDF, which reflects this low temperature 
behavior. As the temperature increases we start seeing a 
third peak in the PDF of the occupation fraction being 
developed close to the ensemble average [see center peak 
in Fig. [T] (b)], so when T — > T g the self-averaging phase 
is approached. 

In Fig. [5] we show the averaged occupation fraction 
versus temperature T using the same deterministic po- 
tential field as in Fig. [TJ For T > T g the average occu- 
pation fraction is Pb(T) and hence as the temperature 
is decreased the average occupation time increases, since 
the particles condensate closer to the minimum of the 
deterministic field which is on x — as the temperature 
is reduced. However when T = T g we see in Fig. [2] a 
type of phase transition in the behavior of the averaged 
occupation fraction, and it does not depend on the tem- 
perature when T < T g , as predicted by our theory Eq. 
(jTTJ) . We note that the convergence of the numerical re- 
sults to the exact ones derived in this paper turned out 
to be relatively slow close to the temperature T = T g if 
compared with temperatures far from T g . 

Finally let us discuss the generality of our results be- 
yond the quenched trap model. We have divided our 
system into two, the observation domain and the rest of 
the system. The partition functions of these domains are 
random variables, due to the randomness of the underly- 
ing Hamiltonian. As we showed, one key ingredient of our 
theory is that the PDFs of these two partition functions 
are one sided Levy stable PDFs. In that case we proved 
that our main Eq. (fl~5| describes the statistics of occu- 
pation time. We know that Levy statistics is common in 
Physics and well established mathematically, and hence 
it seems to us natural to expect that partition functions 
of random systems may have a Levy distribution. In- 
deed a partition function is a sum over energy states and 
if these states are random the connection of the distribu- 
tion of partition functions with limit theorems of sums of 
random variables is expected to be general. In particular 
we could show that our main results describe also occu- 
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FIG. 2: The averaged occupation fraction versus T/T g . 
When T > T g , (p) — Pb(T) namely usual Boltzmann theory 
applies, while for T <T g , (p) — Ps(T fl ) which is independent 
of the temperature T. The lines are theoretical predictions 
Eq. (|17|l and the diamonds are simulation results with no 
fitting. 



pation time statistics in several other models of quenched 
disorder: random comb model which is a simple model 
of a random walk on a loop- less random fractal (l9| . or 
models of anomalous diffusion of a particle on structures 
with distributed dangling bonds in the presence of bias 
2ol 21, 22 1 . Unlike the trap model, in these models ge- 
ometry is the main factor responsible for the non-trivial 
occupation fraction. Sure the exponent does not gener- 
ally turn out to be a — T/T g as in the quenched trap 
model, but our main results are valid as we will discuss 
in a longer publication. 
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